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ABSTRACT 

We report the detection of a magnetic field on the Of?p star HD 108. Spectropolarimetric 
observations conducted in 2007, 2008 and 2009 respectively with NARVAL@TBL and ES- 
PaDOnS@CFHT reveal a clear Zeeman signature in the average Stokes V profile, stable on 
timescales of days to months and slowly increasing in amplitude on timescales of years. We 
speculate that this timescale is the same as that on which Hff emission is varying and is equal 
to the rotation period of the star. The corresponding longitudinal magnetic field, measured 
during each of the three seasons, increases slowly from 100 to 150 G, implying that the polar 
strength of the putatively-dipolar large-scale magnetic field of HD 108 is at least 0.5 kG and 
most likely of the order of 1-2 kG. 

The stellar and wind properties are derived through a quantitative spectroscopic analysis 
with the code CMFGEN. The effective temperature is difficult to constrain because of the 
unusually strong He i A 4471 and He i A 5876 lines. Values in the range 33000-37000 K are 
preferred. A mass loss rate of about 10"^ M© yr~' (with a clumping factor f=0.01) and a wind 
terminal velocity of 2000 km s ' are derived. The wind confinement parameter 77* is larger 
than 100, implying that the wind of HD 108 is magnetically confined. 

Stochastic short-term variability is observed in the wind- sensitive lines but not in the 
photospheric lines, excluding the presence of pulsations. Material infall in the confined wind 
is the most likely origin for lines formed in the inner wind. Wind-clumping also probably 
causes part of the Ha variability. The projected rotational velocity of HD 108 is lower than 
50 km s"', consistent with the spectroscopic and photometric variation timescales of a few 
decades. Overall, HD 108 is very similar to the magnetic O star HD 191612 except for an 
even slower rotation. 

Key words: massive stars - magnetic fields. 



1 INTRODUCTION 

Very little is known about magnetic fields of O stars and very 

few O stars are yet known as magnetic; preliminary results (e.g., 

*„ ^ ^ u ■ . u T.., o JT .-„„Tx iDonati et al.ll2002ll2006l ; rBouret et al.ll2008l : iPetit et alj|2008 f) sug- 
* Based on observations obtained at the Telescope Bernard Lyot (TBL) ' ; " ^' .' 7-" ^ '"- " ^ '' 

J ,,u r^ J c ti •• T r /r^EUTN ,-cxjT ■ , J u gcst that most magnetic O stars could be the high-mass equivalent 

andat the Canada-rrance-Hawaii telescope (CrHf). CrHT IS operated by c o 0-1 

the National Research Council of Canada, the Institut National des Sciences of magnetic chemically peculiar A and B stars (the so-called Ap 

de I'Univers of the Centre National de la Recherche Scientifique of France and Bp stars) with fields of primordial origin, i.e., foss il remnants 

(INSU/CNRS), and the University of Hawaii, while TBL is operated by from the formation stage JDonati & Landstreetl l2009h . They are 

CNRS/INSU. also suspected to be progenitors of highly-magnetic neutron stars 
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and m agnetars (e.g. IPetit et alJIlOOSi : iFerrario & Wickramasinghd 
l2008h . although alternative explanations exist JDavies et alj2009l) . 



Whenever known, rotation rates of magnetic O stars are sig- 
nificantly smaller than non-magnetic O stars of similar spectral 
types, further strengthening the analogy with the lower-mass Ap 
and Bp stars. Evolutionary models of massive stars including mag- 
netic fields also suggest that the internal rotation is strongly affected 
by the presence of magnetic fields, enforcing in particular mostly 
solid-body rotation profiles throughout the bulk of the outer radia- 
tive zone (Maeder & Meynet 2003, 2004). Last but not least, mag- 
netic fields are potentially capable of affecting the powerful radia- 
tive winds of O stars (Babel & Montmerle..l997.) and, in some ex- 
treme cases, of fully dis sipating their angular momentum through 
wind magnetic braking d ud-Doula et al.ll2009l) . Lives of magnetic 
O stars are therefore expected to differ significantly from those of 
their non-magnetic counterparts. 

^In this respect, the prototypical cases of 0' Ori C JDonati et al.l 

l2002h and HD 191612 (Donati et al. 2006) are very interesting. 
The magnetic field of e' Ori C (about 1.1 kG tilted at about 45° 
to the rotation axis) is strong enough to deflect the wind flows from 
both magnetic hemispheres towards the magnetic equator where 
they collide, generating a strong shock, a very hot (coronal-like) 
post-shock region and a denser cooling disc (more or less con- 
fined to the magnetic equatorial plane). Emission in Balmer lines 
exhibits very strong rotational modulation, the amount of emis- 
sion being directly related to the geometry of the cooling disc as 
seen from the Earth: strongest when the disc is seen pole-on an d 
weakest when the disc is seen edge-on (e.g., iDonati et al.ll2002h . 
The presence of magnetic fields and the very similar Ha modula- 
tion reporte d for HD 19 1612 suggests that this star is analogous 
to S' Ori C JDonati et al.l l2006). but more evolved and much more 
slowly rotating (538 d vs. 15.4 d for 0' Ori C) as a likely result 
of wind magnetic braking. Recent (yet unpublished) data (Wade 
et al, in preparation) fully confirms the initial conclusions and the 
remarkable similarity between HD 191612 and 6' Ori C. 

Exploring further the potential impact of magnetic fields on 
the life of massive stars requires new magnetic O stars to be iden- 
tified. From the example of 6/' Ori C and HD 191612, those ex- 
hibiting periodically varying Ho- emission are likely to be very 
good candidates. In particular, HD 108 is interesting. It is a mem- 
ber of the Of?p class composed of the Galactic stars HD 191612, 
HD 148937, NGC 1 624-2 and CPD-28de g 2561, the latter two re- 
cently discovered bv lWalbom et al.l (12010). This class is principally 
defined by emission in the C m 4647-4650 lines at least as strong 
as in the neighbouring N m 4630-4640 lines. Secondary classifica- 
tion criteria include sharp emission or P-Cygni profiles in He i and 
Balmer lines. As 0' Ori C and HD 191612, HD 108 shows spectral 
type variations (from 04 to 08) due to changes in the strength of its 
He I lines. They (and the Balmer lines too) fluctuate between pure 
absorption and P~Cygni profiles. These line morphologies define 
"low" and "high" states respectively. The changes are correlated to 
photometric variations: HD 108 is brighter when the emission com- 
ponents are stronger, again similarly to HD 191612. They seem to 
occur with a period of a few decades JNaze et alj200ll:lBarannikovl 
2007h compared to 5 38 days for HD 191612 JWalbom et aljbooi : 
Howarth et a02007l) . If the analogy with 6)' Ori C and HD 191612 



holds, HD 108 could thus potentially show up as a rather extreme 
case of wind magnetic braking. Note that a marginal detection 
based on a single exposure was reported by iHubrig et al.l 1 120081) 
for the third member of the Galactic Of?p class, HD 148937 (but 
see lSilvester et al.ll2009l , for a discussion of the significance of the 



detection). Further observations are required to confirm this detec- 
tion. 

To investigate this issue, we initiated a spectropolarimet- 
ric monitoring program of HD 108 using NARVAL at the 2 m 
Telescope Bernard Lyot (TBL) in the French Pyrenees and later on 
ESPaDOnS at the 3.6 m Canada-France-Hawaii Telescope (CFHT) 
atop Mauna Kea, Hawaiy. Reports indicate that HD 108 has just 
passed the epoch of minimum Yia emission (Naze et al, in prepara- 
tion), suggesting that longitudinal magnetic fields are also close to 
their minimum value and slowly rising again (provided the analogy 
with 0' Ori C and HD 191612 holds). 

We report here the first results after 3 observing seasons (2007 
to 2009). In particular, we report the detection of a longitudinal 
magnetic field of about 100 G at the surface of HD 108, suggesting 
that the analogy with 6/' Ori C and HD 191612 is likely. We also 
carry out a detailed spectral analysis to obtain updated constraints 
on the parameters of HD 108. We finally discuss briefly the im- 
plications of our results for understanding of how magnetic fields 
affect the lives of massive stars. 



2 OBSERVATIONS 

Spectropolarimetric observations of HD 108 were collected with 
NARVAL® TBL and ESPaDOnS® CFHT in 2007 October, 2008 
October and 2009 July to October. Altogether, 110 sequences (24 
in 2007, 27 in 2008 and 59 in 2009) were obtained in multiple runs, 
each sequence consisting of four individual subexposures taken in 
different polarimeter configurations. 

From each set of four subexposu res we derive a mea n Stokes V 
spectrum following the procedure of Donati et al.l ( 119971) , ensuring 
in particular that all spurious signatures are removed at first order. 
Null polarization spectra (labeled A') are calculated by combining 
the four subexposures in such a way that polarization cancels out, 
allo wing us to check th at no spurious signals are present in the data 
(see Donati et al.ll 19971 for more details on how A^ is d efined). All 
frames were processed using Libre ESpRIT ( IDonati e t al. 1997), a 
fully automatic reduction package installed at TBL and CFHT for 
optimal extraction of NARVAL and ESPaDOnS spectra. The peak 
signal-to-noise ratios per 2.6 km s"' velocity bin range from 140 
to 1 100, depending on the instrument, on the exposure time and on 
weather conditions (see Table[T]for a complete log). 



3 MAGNETIC DETECTION 

Least-Squares Deconvolution (LSD: IDonati et alJll997D was ap- 
plied to all observations. The line list was constructed manually 
to include the few moderate to strong absorption lines that are 
no more than moderately affected by the wind. All lines showing 
strong emission from the wind and/or circumstellar environment at 
the time of our observations (e.g.. Ho-) were excluded from the list. 
The remaining very weak contamination is not circularly polarised. 
Hence it only slightly affects the calculation of the longitudinal field 
but not the Stokes V profile and thus the field detection. The C iv 
lines at 5801.3 and 5812.0 A are used as reference photospheric 
lines from which we obtain the average radial velocity of HD 108 
(about -62 km s"'); about a handful of unblended absorption lines 

' This effort is part of the MiMeS (Magnetism of Massive Stars) interna- 
tional research program aimed at characterizing magnetic fields in upper 
main sequence stars. 
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Table 1. Journal of observations. Columns 1-6 list the date, the instrument used, the range of heliocentric Julian dates, the range of UT times, the number of 
sequences and the exposure time per individual sequence subexposure, and the range of peak signal to noise ratio (per 2.6 km s"' velocity bin), for each night 
of observation. Column 7 lists the rms noise level (relative to the unpolarized continuum level and per 7.2 km s"' velocity bin) in the circular polarization 
profile produced by Least-Squares Deconvolution once averaged over the whole night. Column 8 gives the longitudinal field strength (and the corresponding 
Icr error bar) determined from the night averaged polarization spectrum. 



Date 


Instrument 


HJD 


UT 


'exp 


S/N 


0"LSD 


Bz 






(2,454,000+) 


(h:m:s) 


(s) 




(10-^/,) 


(G) 


2007 Oct 15 


NARVAL 


389.38926-389.47736 


21:14:22-23:21:14 


3 X 4 X 900 


500-590 


2.2 


-94.2±54.3 


2007 Oct 16/17 


NARVAL 


390.48564-390.57375 


23:33:12-01:40:05 


3 X 4 X 900 


350-550 


2.6 


-5L3±63.8 


2007 Oct 18 


NARVAL 


392.33125-392.41936 


19:50:58-21:57:51 


3 X 4 X 900 


340-390 


3.3 


0.9±82.9 


2007 Oct 19 


NARVAL 


393.36308-393.45120 


20:36:51-22:43:45 


3 X 4 X 900 


580-620 


2.1 


-66.6±50.4 


2007 Oct 20 


NARVAL 


394.34587-394.43399 


20:12:07-22:19:01 


3 X 4 X 900 


560-600 


2.1 


-99.6±50.8 


2007 Oct 21 


NARVAL 


395.37385-395.46197 


20:52:28-22:59:21 


3 X 4 X 900 


560-590 


2.1 


-215.1±52.0 


2007 Oct 23 


NARVAL 


397.35103-397.43915 


20:19:43-22:26:37 


3 X 4 X 900 


540-590 


2.1 


-70.6±52.7 


2007 Oct 24 


NARVAL 


398.35262-398.44076 


20:22:03-22:28:59 


3 X 4 X 900 


550-600 


2.2 


-25.6±52.9 


2008 Oct 14 


NARVAL 


754.30940-754.48559 


19:19:21-23:33:04 


5 X 4 X 900 


530-610 


1.6 


-81.0±39.2 


2008 Oct 16 


ESPaDOnS 


755.85707-755.88927 


08:28:03-09:14:25 


2 X 4 X 650 


540-580 


2.7 


-115.4±69.5 


2008 Oct 23 


NARVAL 


763.29085-763.46708 


18:53:06-23:06:53 


5 X 4 X 900 


460-550 


1.9 


-137.0±48.0 


2008 Oct 24 


NARVAL 


764.27121-764.44746 


18:24:52-22:38:41 


5 X 4 X 900 


600-650 


1.5 


-106.9±37.9 


2008 Oct 25 


NARVAL 


765.28297-765.45922 


18:41:52-22:55:41 


5 X 4 X 900 


420-520 


2.0 


-56.4±51.7 


2008 Oct 26 


NARVAL 


766.26993-766.44619 


18:23:09-22:36:59 


5 X 4 X 900 


420-530 


1.9 


-66.5±50.3 


2009 Jul 05 


ESPaDOnS 


1018.04324-1018.10593 


13:00:46-14:31:02 


2 X 4 X 1300 


950-1040 


1.6 


-134.4±43.8 


2009 Jul 07 


ESPaDOnS 


1020.04033-1020.10416 


12:56:22-14:28:16 


2x4x1300 


830-950 


1.7 


-286.0±47.3 


2009 Jul 09 


ESPaDOnS 


1022.04307-1022.10552 


13:00:05-14:30:00 


2x4x 1300 


1070-1090 


1.5 


-222.3±40.1 


2009 Jul 13 


ESPaDOnS 


1026.03827-1026.10162 


12:52:44 14:23:58 


2x4x 1300 


1030-1060 


1.5 


-140.0±40.6 


2009 Jul 14 


ESPaDOnS 


1027.05150-1027.11423 


13:11:41-14:42:01 


2x4x1300 


1100-1110 


1.4 


-142.9±39.8 


2009 Jul 21 


NARVAL 


1033.53871-1033.63538 


00:52:37-03:11:49 


3 X 4 X 750 


330-520 


2.7 


-130.7±66.4 


2009 Jul 25 


NARVAL 


1037.53156-1037.60579 


00:41:56-02:28:49 


3 X 4 X 750 


480-510 


2.5 


-160.0±62.2 


2009 Jul 26 


NARVAL 


1038.57750-1038.64709 


01:47:60-03:28:11 


3 X 4 X 700 


330^50 


3.1 


121.0±82.4 


2009 Jul 27 


NARVAL 


1039.53964-1039.65875 


00:53:23-03:44:53 


4 X 4 X 806 


310^40 


2.7 


-28.1 ±74.0 


2009 Jul 28 


NARVAL 


1040.53824-1040.62633 


00:51:16-02:58:06 


3 X 4 X 900 


260^60 


3.2 


-113.0±87.6 


2009 Jul 29 


NARVAL 


1041.53043-1041.61852 


00:39:56-02:46:46 


3 X 4 X 900 


510-580 


2.2 


-27.8±55.3 


2009 Jul 30 


NARVAL 


1042.54192-1042.63003 


00:56:23-03:03:16 


3 X 4 X 900 


580-640 


2.0 


-62.7±50.7 


2009 Jul 31 


NARVAL 


1043.52399-1043.64394 


00:30:28-03:23:11 


4x4x812 


500-550 


2.0 


-161.6±51.3 


2009 Aug 01 


NARVAL 


1044.54174-1044.62418 


00:55:57-02:54:39 


3x4x812 


540-590 


2.1 


-38.0±55.1 


2009 Aug 02 


NARVAL 


1045.53783-1045.62615 


00:50:14-02:57:24 


3 X 4 X 900 


390^20 


3.1 


-99.4±80.0 


2009 Aug 04 


NARVAL 


1047.56189-1047.64999 


01:24:42-03:31:34 


3 X 4 X 900 


560-600 


2.1 


-161.4±36.6 


2009 Aug 05 


NARVAL 


1048.51890-1048.60701 


00:22:43-02:29:35 


3 X 4 X 900 


580-600 


2.1 


-141.9±52.7 


2009 Sep 03 


ESPaDOnS 


1078.07457 


13:41:03 


4x1300 


940 


2.2 


-109.7±55.8 


2009 Sep 08 


ESPaDOnS 


1082.95060-1083.01300 


10:42:24-12:12:14 


2 X 4 X 1300 


1030-1040 


1.5 


-156.1±38.2 


2009 Sep 25 


ESPaDOnS 


1099.88204-1099.94611 


09:03:29-10:35:43 


2x4x 1300 


990-1010 


1.5 


-208.5±40.3 


2009 Oct 01 


ESPaDOnS 


1105.86247-1105.92525 


08:35:21-10:05:45 


2x4x 1300 


880-960 


1.7 


-140.4±45.2 


2009 Oct 05 


ESPaDOnS 


1109.96249-1110.02567 


10:59:28-12:30:27 


2 X 4 X 1300 


140-600 


3.2 


-273.8±84.1 


2009 Oct 10 


ESPaDOnS 


1114.87331-1114.93701 


08:51:12-10:22:56 


2 X 4 X 1300 


1060-1100 


1.6 


-239.7±41.2 



that are not blueshifted with respect to the reference frame by more 
than 30 km s"' are also included in the list. We ended up with a hst 
of 16 lines, whose characteristics are summarized in Table [2] All 
LSD profiles were produced on a spectral grid with a velocity bin 
of 7.2 km s"', providing reasonable sampling given the significant 
macroturbulent broadening detected on HD 108 (see Sect. 14. 3t . 

Using this line list, we obtained, for all collected spectra, mean 
circular polarization (LSD Stokes V), mean polarization check 
(LSD A') and mean unpolarized (LSD Stokes /) profiles (corre- 
sponding to a weighted - relative to S/N^, line depth and wave- 
length - average line of central wavelength 5500 A and Lande 
factor 1.2). Averaging together all LSD profiles recorded on each 
night of observation (with weights proportional to the inverse vari- 
ance of each profile) yields relative noise levels ranging from 1.4 
to 3.3 (in units of KT'^Ic) in V and A' profiles. In general, no sig- 
nificant magnetic signature was detected in individual spectra. In a 



few exceptional cases, however, nightly averages of spectra showed 
marginal evidence for the presence of magnetic fields at the surface 
of HD 108 (see column 8 of Table[T](. 

Using a line list containing the five lines with weakest inter- 
stellar contamination and smallest blueshifts with respect to the 
stellar rest velocity (indicated in bold face in Table|2ll, the detected 
Zeeman signature is compatible (within noise level) with that de- 
rived using the full line list, confirming that the detected signature 
is due to magnetic fields at the surface of the star. 

The detected Zeeman signature (featuring a peak-to-peak am- 
plitude of 0.04%) cannot be attributed to the small crosstalk be- 
tween circular and linear polarisations that ESPaDOnS and NAR- 
VAL used to sufi'er (smaller than 1% on ESPaDOnS during our 
2009 run, smaller than 4% for NARVAL and for ESPaDOnS at 
all other epochs) as this would imply that photospheric lines of 
HD 108 exhibit roughly antisymmetric linear polarisation signa- 



4 Martins et al. 



Table 2. Lines used for Least-Squares Deconvolution. The line depths (col- 
umn 3) were directly measured from our spectra while the Lande factors 
(column 4) were derived assuming LS coupling. Bold face indicates lines 
with the weakest wind contamination 



Wavelength 


Element 


Depth 


Lande 


(A) 




(/c) 


factor 


4026.198 


He I 


0.33 


1.1 


4199.839 


Hen 


0.15 


1.0 


4379.201 


Nni 


0.06 


1.1 


4387.929 


He I 


0.08 


1.0 


4471.483 


He I 


0.41 


1.1 


4510.963 


Nra 


0.06 


1.1 


4514.854 


Nni 


0.08 


1.2 


4541.593 


Hen 


0.20 


1.0 


4713.139 


He I 


0.13 


1.5 


4921.931 


He I 


0.14 


1.0 


5015.678 


He I 


0.11 


1.0 


5411.516 


Hen 


0.25 


1.0 


5592.252 


Oiii 


0.40 


1.0 


5801.313 


Civ 


0.20 


1.2 


5811.970 


Civ 


0.15 


1.3 


7065.176 


He I 


0.22 


1.5 



HD 108, LSD profiles, 2009 
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Figure 1. LSD circularly-polarised (Stokes V), polarisation check (N) and 
unpolarised (Stokes /) profiles of HD 108 (top/red, middle/green, bot- 
tom/blue curves respectively) averaged over all NARVAL and ESPaDOnS 
data collected in 2009. LSD V and A' profiles are expanded by a factor of 
250 and shifted upwards (by 1.15 and 1.05) for display purposes. 



tures with peak-to-peak; amplitudes of at least 1-4% - a highly un- 
likely possibility. It would also imply that the Zeeman signature 
detected with ESPaDOnS in 2009 is smaller (by at least a factor of 
4) than those measured with NARVAL and ESPaDOnS at all other 
epochs, reflecting the instrumental fix implemented on ESPaDOnS 
in early 2009; since this is not the case, we can safely claim that the 
detected Zeeman signature is not caused by a crosstalk problem. 

For each of the main observing seasons (2007, 2008 and 
2009), we constructed average LSD signatures and corresponding 
variances. In 2009, the mean Stokes V signature is clearly detected 
with a reduced chi-squared x] within the line profile of about 12|j 



^ The reduced chi-squared estimator quantifies how significantly the ob- 
served Zeeman signature departs from a null profile (i.e. V=0).) 



(see Fig. [TJ; in 2008 and 2007, the corresponding xl values (2.5 
and 1.5 respectively) are smaller but large enough to claim a defi- 
nite detection in 2008 and a marginal one in 2007 (with a 98% con- 
fidence level). No signal is detected in the A' (polarization check) 
profile, demonstrating that the signal we detect is truly attributable 
to circular polarization. The variance profiles show no evidence for 
variability of the Zeeman signatures throughout each season, even 
in 2009 where the signal to noise ratio is highest. We therefore 
conclude that the Zeeman signature, if varying, is changing very 
slowly on typical timescales of years rather than months or weeks. 
The longitudinal fields we measure in 2007, 2008 and 2009 (us- 
ing the first order moment method of lDonati et alJI 19971 . - inspired 
by Semel 1967 - and integrating through a velocity range of -300 
to 120 km s"') are respectively equal to -90 ± 20, -100 ± 20 and 
-150± 10 G (Icr error bars), suggesting a slow rise in field strength 
on timescales of years. 

The detection of a magnetic field on HD 108 as well as the 
lower limit of months to years on the timescale on which the Zee- 
man signatures are varying both strengthen the aforementioned 
analogy with 6' Ori C and HD 191612. For those stars, the lon- 
gitudinal fields and line emissions (and photometric brightness 
for HD 191612) are all modulated with the same period (the ro- 
tation period). We thus speculate that the longitudinal field of 
HD 108 tightly correlates with line emission and photospheric 
brightness, weakest at minimum emission and minimum brightness 
and strongest at maximum emission and maximum brightness (as 
in HD 191612). This is in agreement with the slow increase in both 
longitudinal field and Balmer line emission that we observe from 
2007 to 2009. 

If we assume that the magnetic field at the surface of HD 108 
is grossly dipolar (e.g., like that of 0' Ori C and HD 191612, 
iDonati et al. 2002. .2006), it implies that the polar field at the sur- 
face of the star (r oughly 4 times larg er than the maximum longi- 
tudinal field, e.g.. lDonatietalJl2002l) is at least 0.5-0.6 kG. This 
estimate is bas ed on the clas sical relati on between longitudinal an d 
polar field of lPrestod ^1967* ) (see also lDonati & Landstree^l2009h . 
If we further take into account that HD 108 is still at least 10 yr 
away from longitudin al field maximum (e.g. given the observed 
photometric variations, lBaraiuiikovl2007l) , we conclude that the po- 
lar field of HD 108 is likely stronger than 1 kG and potentially as 
strong as (or even stronger than) that of HD 191612. 



4 STELLAR AND WIND PHYSICAL PROPERTIES 



iNaze et al.l J2008I) conducted a study of the "high" optical state of 
HD 108. Here, we present a new analysis based on our NARVAL 
spectra collected in 2008 and 2009, as well as FUSE data obtained 
in 2007 and 2008 and retrieved from the MAST archive. These data 
correspond to the "low" state of HD108. The spec tral analysis was 
perfo rmed using the atmosphere code CMFGEN JHillier & Milleil 
Il998h . 

An exhaustive description of the atmosphere models can be 
found in lHillier & Milleii jl998r) . In practice, an iterative scheme is 
developed. At each step, the radiative transfer equation is solved 
in spherical geometry. With the resulting specific intensity, the 
rate equations yield non~LTE level populations which in turn are 
used to compute opacities. These opacities are subsequently in- 
jected in the radiative transfer equation in the next iteration step. 
The temperature structure results from the radiative equilibrium 
equation. Line-blanketing is included by means of the super- 
level approach. The atmosphere velocity structure is given as in- 
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put and results from the connection o f a pseudo-hydrostatic struc- 
ture (taken from TLUSTY models, iLanz & Hubeiivl l2003h with 
a so-called /3-velocity law representing the wind structure (e.g. 
iLamers & Cassin elli 1999). The density structure directly follows 
from mass conservation. All models take into account clump- 
ing, a nd Auger ionizat ion by X-rays (using the X-ray luminos- 
ity of Naze et al. 20041) . We adopt a solar chemical composition 
JGrevesse et alj|2007 ) except for helium and nitrogen (see below). 
Once the model atmosphere is obtained, a formal solution of the 
radiative transfer equation is performed with to produce the final 
synthetic spectra. 

The derived stellar and wind proper t ies of HD 108 are gath- 
ered in Table [3] The results of lNaze et alj ( 120081) are also given for 
comparison. 



4.1 Photospheric parameters 

The effective temperature Teff is usually derived from the relative 
strength of He i to He n lines in O stars. In the case of HD 108 
however, this procedure turned out to be considerably more com- 
plex. First, some He i lines present a central emission component 
the origin of which is not known (e.g. He i A 4144, He i A 4388, and 
He I A 4920). Second, most lines have highly asymmetric profiles, 
their blue wing rising towards the continuum less steeply than their 
red wings (see e.g. Fig.[T](. Third, He i A 5876 and, more worrisome, 
the normally well-behaved He i A AAl 1 line show unusually strong 
absorption in the 2007-2009 spectra. In fact, these two lines are the 
strongest features of the entire optical spectrum. Fig. [2] illustrates 
the problems encountered when fitting the He i lines. It shows the 
difference between the observed He i equivalent widths and those 
of various models. Clearly, He i A AA1\ and He i A 5876 are always 
too weak in our models, because both the line depth and equivalent 
width are weaker than observed. The situation is improved only for 
very high He/H ratios (> few tens), but such values are unrealistic 
and provide a worse fit to other He lines. We thus decided to leave 
these lines aside and focused on the remaining He i and He n lines. 
The best fit was achieved for T^ff = 33000-37000 K and a He/H 
ratio of 0.3 (by number). It is shown in Fig. [2] with Teff = 35000 
K. We also needed a larger than solar N abundance to fit the N m 
lines in the interval 4510-4540A, namely, N/H = 3.6 x 10"* (~6x 
solar by number). A discussion of the abundance pattern of mag- 
netic O stars in relation with evolutionary models will be presented 
in a subsequent publication (Escolano et al. 2010, in prep). 

HD 108 is a member of the OB association Cas 0B5. The lu- 
mino sity L^ was determ ined by assuming a distance of 2.51 ±0.15 
kpc JHumphrevsl 119781) and fitting the spectral energy distribu- 
tion using lUE data and UBVJHK photometry. The uncertainty on 
log j^ is of about 0. 1 dex and is entirely dominated by the uncer- 
tainty on the distance. The surface gravity (log g) was obtained 
from the fit to the Wfi and Hy line wings. A value of 3.5 (typical 
of late type giants/supergiants; see lMartins et alj|2005l) is derived, 
with an uncertainty of about 0.2 dex. The position of HD 108 in the 
HR diagram is shown in Fig. [3] HD 191612 and 0' Ori C are also 
included for comparison purpose. HD 108 is about 4 Myr old just 
as HD 191612, reinforcing the analogy between both stars. Given 
its slightly higer mass, HD 108 might be a little more evolved than 
HD 191612. HD 108 is also older than e' Ori C, consistent with 
the scenario according to which 0' Ori C represents a precursor of 
HD 191612 toonati et al J2006i) and thus HD 108. 

The comparison of the results of the present study to the one 
bv lNaze et alj yOOSJ) shows that the derived stellar parameters are 



Table 3. Summary of stellar and wind properties of HD 108. 

present work Naze et al. 2008 



T.//(K) 


35 000 ± 2000 


37000 ± 2000 


log g (cgs) 


3.50 ± 0.20 


3.75 ±0.10 


log L*/Lo 


5.70 ±0.10 


5.40 ±0.10 


R*/Ro 


19-2 !^i 


12.3 ±2.1 


M*/Mo 


43 +32 


-35 


\sini (km s"') 


<50 


-40 


logM 


.7 O+O-O 
'■"-0.40 


-6.5/-7.0 


Vec (kms"') 


2000 ± 300 


2000 


Vc/ (clumping; km s"') 


30 


30 


/ (clumping) 


0.01 


0.01 


E(B-V) 


0.47 


- 


log LxILbol 


-6.2 


- 


He/H (by number) 


3.0x10-' 


o.r 


N/H (by number) 


3.6x10-* 


- 



in rather good agreement in spite of the difficulties encountered in 
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the present analysis. Only the luminosity is marginally higher in the 
present study (b ut the method we u sed for its determination is more 
accurate than in lNaze et alj|2008h . This overall similarity demon- 
strates that the variability altering the shape of the optical lines is 
probably not due to changes in the star's stellar parameters. Instead, 
they might be related to the wind or immediate environment of the 
star. 



4.2 Wind parameters 

We relied on strong EUV and FUV lines to constrain the mass- 
loss rate (M) and the terminal velocity (Vo<,). We avoided the use of 
wind sensitive optical lines (Hce and He ii A 4686) since they show 
a high degree of short-term variability and can be contaminated by 
non-stellar emission. 

We used only the FUSE LiF2A channel of the FUSE spectra 
in our analysis since wavelengths shorter than llOOA are heavily 
contaminated by the ISM. Our best fit is shown in Fig. |4] (upper 
panel) and the corresponding parameters are listed in Table [3] The 
synthetic P \ AA 1118,28 and Si iv A/ll 122,28 lines are somewhat 
weaker and narrower than observed. By decreasing the mass-loss 
rate we tend to decrease this discrepancy, but a worse fit is achieved 
to the C III /11176 and C iv /ill69 transitions. 

Unfortunately, no FUV spectra (1200 to 2000 A) exist for the 
"low" optical state. The only lUE observations cover the period 
1978-1989 and are quite similar, showing no signs of significant 
variability. They presumably correspond to HD 108's "high" op- 
tical state (the Bal mer lines being i n emission at lea st between 
1987 and 1989, see IUnderhilllll994 iNaze et alj|200l[) . However, 
we found out that the set of parameters derived from the FUSE 
spectra lead to a good fit of the lUE data (see lower panel of Fig. 
O. We thus conclude that the global wind properties of HD 108 do 
not vary significantly between the low and high state. In particu- 
lar, the mass loss rate does not differ by more than a factor of two 
(corresponding to the uncertainty on our determination). We note 
that the wind par ameters are also consistent with those derived by 
iNaze et al.l ( 120081) who used the same lUE data but different optical 
spectra. 

Given that the wind properties do not seem to depend on the 
star's state (high or low), we used N iv /1 17 18 to constrain th e 
clumping factor (see a complete description in lBouret et alj|2005n . 
A small / value (0.01) corresponding to a large degree of wind in- 
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Figure 2. Determination of the effective temperature and surface gravity of 
HD 108. Upper panels: final model from CMFGEN (red/grey) and observed 
spectrum (black). Lower panel: observed minus synthetic equivalent widths 
of several Helium lines (AWa (O-M) = Wf ' - W™"''''). A perfect agreement 
is illustrated by the dash-dotted line. 



homogeneities was necessary to correctly reproduce the blue wing 
extension. 



4.3 Rotational velocity and macroturbulence 

We used the Fouri er transform method to constrain the projected 
rotational velocity ( lGravlll992h . To avoid any problem caused by 
wind and/or circumstellar emission, we selected photospheric lines 
showing the least degree of contamination. It turned out that no line 
showed a pure symmetric profile. Even for the best candidate fea- 
tures, the blue wing clearly showed a wider extension than the red 
wing, indicating some kind of contamination, most likely due to 
the wind. Consequently, we decided to work on an artificial, sym- 
metric line profile created from the red wing of C iv /I 5812, which 
we assume better reflects the true photospheric properties. We used 
the average of the July-August 2009 NARVAL spectra as input. No 
minimum was seen in the Fourier transform above the noise level, 
indicating that 1) the maximum vsini compatible with our artificial 
profile is ~ 50 km s"' and 2) that macroturbulence is an impor- 
tant broadening mechanism. Tests run on N iv /t 4057 and O m A 
5592 gave similar results. Note that the origin of macroturbulence 
in O stars is widely unknown, although Aerts et al. (2009) recently 
claimed that non-radial pulsations could be an important ingredi- 
ent. 

To further investigate the rotational velocity and to quantify 
the amount of macroturbulence, we used a synthetic line profile 
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Figure 3. HR diagram with the position of HD 108 (red circle), HD 191612 
(blue triangle) and 0* Ori C (green asterisk) indicated. Evolutionary tracks 
include stellar rotation and are from Mevnet & Maeder (2003). Isochrones 
are sho wn as dotted lines. Ste llar parameters for HD 191612 and S' Ori C 
are from lWalbom et al.l J2003h and lSimon-Diaz et al.l 12006) respectively. 
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1300 1400 1500 1600 1700 
Wavelength (A) 

Figure 4. UV and far-UV observed spectra (black) and our final CMFGEN 
model (red). Top: FUSE LiF2A spectrum. Bottom: WE spectrum. The re- 
gion below 1200 A common to the FUSE spectrum is very noisy and is not 
shown. The flux is in units of erg s"' cm"^ A"' . 
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(taken from the TLUSTY OSTAR2002 grid of models). We arti- 
ficially scaled it in order to match the equivalent width of the ob- 
served C IV /i 5812 line. Since macroturbulence is the main broad- 
ening mechanism (see below) this is justified, the intrinsic profile 
being unimportant. We subsequently convolved this synthetic pro- 
file with 1) a rotational profile (characterized by vshi;) and 2) a 
pure Gaussian profile (characterized by its FWHM and the corre- 
sponding macroturbulence velocity Wrn-ni)- The latter mimics the ef- 
fect of isotropic turbulence. Several (vszn/ / FWHM) combinations 
were used and compared to the artificial profile described above. 
We found that rotational velocities larger than about 50 km s"' 
are indeed excluded and that a macroturbulence larger than ~ 30 
km s"' was necessary to reproduce the wings slope and extension. 
Different combinations of {ysini / Vmac) gave fits of similar qual- 
ity. Fig.|5]shows three examples with v&ini - 0, 50 and 80 km s"' 
and the corresponding Vj^ac = 45 , 38, km s~'. HD 108 is there- 
fore very similar to HD191612 JHowarth et alj|2007l) which also 
features macroturbulence-dominated line profiles. 



5 DISCUSSION 

The analogy between HD 108 and (9' Ori C/HD 191612 suspected 
in Sect.[T]seems to be confirmed by the present study. In particular, 
the similarity with HD 191612 is striking. Both stars show the same 
long-term spectroscopic variability, with He i and Balmer lines 
(other than Ha) chang ing from P-Cygn i to absorption-domina ted 
profiles (see Table 6 of lNaze et alj200lh . lHowarth et al.l ( l2007h re- 
ported a strong correlation between Ha emissi on and photometry , 
HD 191612 being brighter in the "high" state. iBarannikovl IJOOTh 
provides evidence for the same behaviour of HD 108: photometry 
over the last 15 years reveals a fading of HD 108 in parallel to a 
decrease of the emission component of He i and Balmer lines (see 
also Naze et al. in prep.). 

A notable difference is the variability timescale. While a pe- 
riod of 538 days is identified for HD 191612, a value of several 
decades is preferred for HD 108. Thi s rough estimate re lies on two 
facts. First, the photometric survey of lBarannikovl ( 120071) discussed 
above points to a timescale of at least 15 years (and of at least 
30 years assuming roughly sinusoidal photometric variations). Sec- 
ond, a transition between absorption/P-Cygni of several line pro- 
files was observed at least twice during the last c entury, pointing t o 
a periodicity of abo ut 50-60 year s (see T able 6 of lNaze et al.l200ll) . 
If, as suggested bv lDonati et"al] (12006^ for HD 191612, the vari- 
ability timescale corresponds to the rotation period, then HD 108 
would be an extremely slow rotator with an equatorial velocity <0. 1 
km s~' (much lower than pred icted by evolutionary models, see Fig. 
1 of lMevnet & Maedenl2003h . Our upper limit on the projected ro- 
tational velocity is fully consistent with that scenario. One might 
thus suspect that HD 108 experienced a strong magnetic braking 
during its evolution. 

Using the stellar and wind parameters derived in Sect. |4] and 
the minimum polar field discussed in Sect. [3] we obtain a con- 
finement parameter rj-,, = j^ > 100. Such a value im plies that 
the w ind of HD 108 is magnetically confined Cud-Doula & Owockil 
l2002t) . If the polar field is as large as that of HD 191612 (1.5 kG, 
see Sect. [3}, this would imply r]^, ~ 800 and thus a strong wind 
confinement. 

According to lud-Doula et alj ( 120091) (see also discussion by 
iDonati et alj2006l) the typical spin-down timescale due to magnetic 
braking of a hot star with a strongly magnetically confined wind is 
Tspin = ^^~R= (k is a constant with a typical value of 0.1 to 0.3). 



For a kG polar field this means Tj,,,,, ~ 1 - 3Myr. Given its Teff and 
luminosity, HD 108 is 3 to 5 Myr old (see Fig.[3l(. According to 
this simple estimate the rotation rate of HD 108 could have been 
reduced by a factor 3 to 200 just by magnetic braking (implying an 
initial rotational velocity of a few to 20 km s"'). However, tailored 
simulations relying on additional observations providing a better 
characterization of the field strength and geometry are obviously 
required to quantitatively tackle the q uestion o f magnetic braking 
iMaeder & Mevnet 2005 : ud-Doula et alj2009l) . 

The spectrum of HD 108 is varying on timescales of a few 
decades. Our excellent time coverage in 2009 allows us to investi- 
gate the existence of faster variability. Fig. |6] shows a selection of 
He I, He ii and Balmer lines t ogether with the cor responding Time 
Variance Spectrum (TVS, see lFullerton et alj|l996.) which quanti- 
fies the degree of variability. All He i and Balmer lines are variable 
on timescales of days. He ii A 4686 also shows fluctuations. But the 
He II A 4200 and He ii A 4542 lines (as well as C iv AA 5801-5812 
and O III/15592, not shown in Fig.|6j are stable over time. The He i 
and Balmer lines, as well as the wind line He ii A 4686, are formed 
over a large fraction of the atmosphere (i.e. photosphere -I- wind), 
especially if a disk/equatorial density enhancement exists. On the 
contrary. He ii A 4200 and He ii A 4542 are formed in the hottest 
part of the atmosphere, in or very close to the photosphere. The 
observed variability thus originates in the wind. Non-radial pulsa- 
tions can thus be excluded since they would trigger variability in 
photospheric lines. A periodicity search performed on equivalent 
width measurements of several lines and using the Lomb-Scargle 
formalism did not reveal any pattern. We thus conclude that HD 108 
presents day-to-day stochastic wind variability. 

In the standard magnetically confined wind scenario, and in 
the case of moderate to strong confi nement, a disk is formed i n 
the inner atmosphere (e.g. Fig. 3 of ud-Doula & Owocki 120021) . 
Thi s disk is denser in its iimer parts. lud-Doula & Owockil ( l2002h 
and ud-Do ula et alj ( 120060 have shown that in absence of sufficient 
centrifugal support, material accumulated in the disk and located 
below t he co-rotation radius falls b ack onto the stellar surface (see 
Fig. 4 of lud-Doula & Owockil2002l) . This process is stochastic, and 
the disc in maintained dynamically by new material injected via the 
channelled wind. In most lines of Fig.|6|we see that variability is lo- 
cated preferentially on the red part of the profile. This can be inter- 
preted as a direct evidence for the infall of "blobs" of disk material. 
Indeed, density fluctuation along the line of sight will cause varia- 
tion of the line strength, and the receding velocity will shift these 
variations towards longer wavelengths. We think that this variabil- 
ity in the red wings of most line profiles is a n indication of mat erial 
infall. Similar conclusions were reached by I Wade et al.l ( 120060 for 
S' Ori C. Ha seems to deviate from this scenario: its variability is 
located over the entire profile. Contrary to most lines of Fig. [6] 
Ha is formed further away from the photosphere and is usually 
more affected by the wind. It is known to be variable in a number 
of non-magnetic O stars, mainly due to the presence of inhomo- 
geneities ("wind-clumping") most likely due to hydrodynamical 
instabilities (e.g. Runacres & Owocki 2002). Markova et al. (200^ 
showed that such clumps could explain the Ha variability observed 
over the entire line profile of O supergiants. We thus suggest that 
in HD 108, Ha is formed in a zone where the putative disk of con- 
fined material is more tenuous than in the formation region of all 
other lines of Fig. [6] In that region, wind-clumping affects the line 
profile in addition to material infall. This explains that variability 
is not confined to the red wing but is located over the entire profile. 
The large clumping factor derived in Sect. |4.2| (f=0.01) is consistent 
with a very inhomogeneous wind. Obviously, more information on 
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Figure 5. Comparison of the C iv /I 5812 line with synthetic profiles with the following vsmi'/Vmac values (in km s"'): 80/0 (left), 50,38 (middle), 0/45 (right). 
The black solid line is the observed profile and the red dashed line the synthetic profile. While vsini larger than 50 km s"' is excluded, any value below this 
limit gives a satisfactory fit to the red wing provided a significant amount of macroturbulence is included. 



the field strength and geometry are needed to better constrain the 
wind confinement and to test our suggestion. In particular, a better 
knowledge of the position of the Alfven radius relative to the Ha 
formation region is important to test our hypothesis that Ha vari- 
ability is at least partly caused by wind-clum ping. 

It is al so worth noting that th e models of lud-Doula & Owockil 

( 120021) and lud-Doula et alj ( 120061) consider only the case of mag- 
netic axis aligned with the rotation axis. In our case a tilt is likely 
present given the long-term line variability, so that the spin-down 
timescale and centrifugal support cannot be strictly compared to 
the theoretical predictions. 

We finally stress that in the case of magnetically confined 
winds, one might question the use of ID atmosphere models to 
constrain the wind properties. However, it is unlikely that our mass 
loss rate estimate is wrong by orders of magnitudes. Indeed, even if 
half of the material initially blown by radiative acceleration was to 
fall back onto the stellar surface, strong signatures of infall such as 
inverse P-Cygni profiles or redshifted absorption should be seen. 
Since the only evidence for infall is the line variability of absorp- 
tion lines, we think our estimate of the mass loss rate is realistic. 



6 SUMMARY AND CONCLUSIONS 

We have presented spectropolarimetric observations of the Of?p 
star HD 108 conducted with the NARVAL and ESPaDOnS in- 
struments (at TBL and CFHT respectively). One hundred and ten 
circularly-polarized spectral sequences have been collected be- 
tween 2007 and 2009. We report the clear detection of a Stokes 
V Zeeman signature stable on timescales of days to months, but 
likely slowly increasing in amplitude on timescales of years. We 
speculate that this timescale is the same as that on which the Ha 
emission is varying and is equal to the rotation period of the star. 
The corresponding longitudinal magnetic field is of order of 100- 
150 G, implying that the polar strength of the putatively-dipolar 
large-scale magnetic field of HD 108 is at least 0.5 kG and most 
likely of order of 1-2 kG. 

The stellar and wind parameters of HD 108 have been derived 
through atmosphere modelling with the code CMFGEN. The ef- 
fective temperature was poorly constrained due to the surprisingly 
strong He i A 4471 and He i A 5876 lines. Values of Teff in the range 



33000-37000 K were conservatively derived from the analysis of 
the remaining He i and the He ii lines. A mass loss rate of the or- 
der 10"' Mq yr"' was derived. The wind is found to be strongly 
clumped (f=0.01). It is also significantly confined (rj-,, > 100, pos- 
sibly up to 800) by the magnetic field. HD 108 shows long-term 
variability in most He i and Balmer lines, with profiles changing 
from pure absorption to P-Cygni on a few decades period. We sug- 
gest that this periodicity corresponds to rotational modulation. The 
(magnetic) equatorial density enhancement implied by the wind 
confinement would then be seen edge-on at minimum line emis- 
sion, and pole-on at maximum emission. The timescale for mag- 
netic braking computed with the derived stellar, wind and magnetic 
properties is consistent with the star's age. A slow rotation is also 
consistent with the low vsini we derive. HD 108 might thus be an 
even more extreme case of slowly rotating magnetic O star than 
HD 191612. 

Short-term line profile variability is also observed in He i and 
wind sensitive lines. Photospheric He ii lines are very stable. Lines 
formed just above the photosphere show evidence for infall while 
the variability of lines formed in the outer wind is more typical of 
normal hot stars' outflows. We suggest that in the inner wind, we 
are witnessing the infall of material channeled by the field lines to 
the magnetic equator and subsequently pulled back to the stellar 
surface by gravity. 

Most of the suggestions discussed in Sect.|5]and summarized 
above are still speculative. Future monitoring of HD 108 is obvi- 
ously needed to confirm the expected coiTelation between the vari- 
ation of the longitudinal magnetic field and the long-term spec- 
troscopic variability. If the few decades spectroscopic and pho- 
tometric modulation corresponds to the rotational period, a com- 
plete mapping of the magnetic topology is not possible until a few 
decades. However, crucial information regarding the field strength 
and geometry can be gathered in the next years since we have just 
passed the phase of minimum emission. The field is thus expected 
to strengthen, making detection easier. Further constraints on the 
field morphology will help to see 1) if the slow rotation of HD 108 
is really a due to magnetic braking and 2) if the short-term variabil- 
ity is caused by dynamical phenomena predicted by current simu- 
lations of magnetically-channelled winds. 



Detection of a magnetic field on HD 108 9 






-200 200 

velocity [km/sj 



(a) Ha 



(b)H8 



(c)H7 






(d) He 1^4471 



(e) He 1/14713 



(f) He I i 4920 




^ 0.004 

: 0,002 




}Vtf*(t#i^,;VV*Mf^^ 



^ 0.004 [r 
I 0.002 



ii^^it^jj^fiiiijirn:^^ 




(g) He II A 4200 



(h) He II A 4542 



(i) He II A 4686 



Figure 6. Line profile variability of Ha, Ufl, Hy, Hei /i 4471, Hei A 4713, Hei A 4920, Hen A 4200, Hen i 4542 and Hen A 4686 (panels a to i) between July 4"' 
and October 9''' 2009. The red bold line is the 2009 average spectrum while the thin lines correspond to individual spectra. All spectra are shown in velocity 
space around the respective rest wavelength. The temporal variance spectrum (TVS, see FuUerton et al. 1996) is also given for each line (lower panels). The 
dot-dashed fine shows the one tr deviation of the continuum, while the dotted fines shows the 99% confidence limit. Variabifity is clearly detected in all Balmer 
and Hei lines. On the contrary, the Hen lines are stable, with only a small variability in the wind line Hen A 4686. 
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